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Thermal activation and macroscopic quantum tunneling in current-biased discrete Josephson 
transmission lines are studied theoretically. The degrees of freedom under consideration are the 
phases across the junctions which are coupled to each other via the inductances of the system. The 
resistively shunted junctions that we investigate constitute a system of A'^ interacting degrees of 
freedom with an overdamped dynamics. We calculate the decay rate within exponential accuracy 
as a function of temperature and current. Slightly below the critical current, the decay from the 
metastable state occurs via a unique [rigid) saddle-point solution of the Euclidean action describing 
the simultaneous decay of the phases in all the junctions. When the current is reduced, a crossover 
to a regime takes place, where the decay occurs via an elastic saddle-point solution and the phases 
across the junctions leave the metastable state one after another. This leads to an increased decay 
rate compared with the rigid case both in the thermal and the quantum regime. The rigid-to-elastic 
crossover can be sharp or smooth analogous to first- or second-order phase transitions, respectively. 
The various regimes are summarized in a current-temperature decay diagram. 



I. INTRODUCTION 

The decay of systems in a metastabla-state has been 
intensively studied in the last decades.EJ At high tem- 
peratures the decay is induced by thermal activation. 
Lowering the temperature, under well-defined conditions 
quantum fluctuations become important and a crossoiwei; 
from thermal activation to quantum tunneling occurs .Bu 
For systems with several interacting macroscopic degrees 
of freedom (DOF), the decay scenario is quite complex. 
For example, a crossover from rigid decay, where all DOF 
decay instantly, torelastic decay, where the DOF decay 
one after anotheioia xan take place in systems with two 
DOF i^ii^he thermalQ as well as in the quantum decay 
rcgimcDEl Recently, the rigid-to-elastic crossover was in- 
vestigated theoreticallyi-jn systems with N DOF in the 
high-temperature limit B In this work we study the var- 
ious decay regimes that occur in these systems at low 
temperatures, where quantum tunneling is relevant. 

A system in a metastable state decays most probably 
via the lowest-lying saddle-point configuration of the Eu- 
clidean action. A crossover in the decay rate occurs, if 
a new lower-lying saddle point appears upon tuning an 
external parameter. It can be continuous or discontinu- 
ous. A continuous or second-order crossover occurs when 
a saddle-point bifurcation takes place. For example, at 
the crossover from thermal to quantum decay, the static 
saddle point that determines the decay at high tempera- 
tures bifurcates into dynamic extrema (instantons) that 
have a lower Euclidean action. This causes an enhance- 
ment of the escape rate F for temperatures T below the 
crossover temperature Tq (see Refs. [2,3]). In steepest- 



descent approximation r(T) and its derivative r'(r) are 
continuous at Tq, whereas the second derivative T"{T) di- 
verges. In analogy to a classical phase transition, Larkin 
and Ovchinnikov denoted it a second-order crossover.0 
They pointed out that for some potential shapes Y'{T) 
can become discontinuous at Tq and the crossover is then 
of first order. Transitions of this type were later investi- 
gated in more detail by ChudnovskiJl3 and studied in var- 
ious physical systems by others.tllE3 The crossover from 
rigid-to-elastic decay is similar and both first-order and 
second-ordsi crossovers have been found in systems with 
two DOFEB 

The DOF we study in this paper are macroscopic 
and coupled to their environment. At low tempera- 
tures, when macroscopic quantum tunneling (MQT) oc- 
curs, the interaction of the DQF with the environment 
leads to quantum dissipation.Ej Among the experiments 
that have been successfully interpreted in the frame- 
work of MQT with dissipation are the investigations of 
the decay of metastable states in current-biased Joseph- 
son junctions (JJ's) and superconducting quantum in- 
terference devices (SQUID's)£il The relevant collective 
coordinate in the case of the JJ is the phase difference 
across the junction. The dc SQUID consists of two par- 
allelly connected JJ's and thus enables the study of the 
quantum dynamics of two coupled degrees of freedom, 
the phases across the junctions. The natural general- 
ization to N degrees of freedom are parallel coupled 
one-dimensional Josephson-junction arrays, also known 
as discrete Josephson transmission lines (DJTL's). The 
continuum limit N ^ oo corresponds to a long JJ if 
the length of the junction I is larger than the Josephson 
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length Aj. The thermaL|decay of the phase in under- 
damped long JJ's (Refs.ErtZl) and overdamped DJTL's 
(Reffl) has been investigated recently. The quantum tun- 
neling in underdamped long JJ's has been analyzed in 
the limit I oo (Ref. [15]) and in the rigid regime l^lc, 
where Ic ^ Aj is the critical length ahave which elastic 
decay via boundary nucleation sets in.llll To date, the in- 
teresting case Z>/c has not been treated in the quantum- 
tunneling regime, neither in the overdamped nor in the 
underdamped case. Solving the problem is difficult since 
then the saddle-point solutions of the action are inho- 
mogeneous in time and space. Two reasons motivated 
us to study N harmonically coupled DOF trapped in a 
metastable state. First, our model can be used to de- 
scribe the decay of the JJ phases in a DJTL and in the 
continuum limit — > oo to represent a long J J. The sec- 
ond reason is the need to analyze the so far undiscussed 
overdamped limit. 

Since the decay close to the rigid-to-elastic crossover 
is determined by the long-wavelength excitations of the 
system, much insight can be gained from studying jdr; 
SQUID'S. Experiments on dc SQUID's (Refs.ErEJ) 
showed that at high currents the classical decay mainly 
proceeds via the simultaneous activation of both phases 
across a common saddle point in the potential-energy 
landscape. However, below a crossover current, it was 
found that this sa|ddle point splits and the phases decay 
one after another A similar phenomenon was proposed 
to occur in the quantum decay process in underdampedQ 
and overdampedo SQUID's, known as the instanton split- 
ting. In this case, the common saddle point of the Eu- 
clidean action bifurcates into two saddles of lower action. 
In addition, a regime was found, where upon tuning tem- 
perature or current a first-order transition takes place 
before the saddle-point bifurcation occurs. 

The rigid-to-elastic crossover in DJTL's is similar. Of 
course, instead of two DOF, one then has A^ coupled 
DOF and the theory has to be generalized accordingly. 
A perturbative treatment to calculate the elastic instan- 
ton solutions of the Euclidean action has been sketched in 
Refs. [7,8] for case of the dc SQUID. In the present publi- 
cation, the perturbation scheme is systematically worked 
out and applied to the problem of an overdamped DJTL 
in a metastable state. With this procedure one is able 
to calculate the split-instanton solution to arbitrary pre- 
cision. We present quantitative results for the quantum 
decay rate of overdamped dc SQUID's, DJTL's and long 
JJ's close to the the rigid-to-elastic crossover. We further 
construct the decay diagram of DJTL's and long JJ's and 
compare it to the one of dc SQUID's. 

The paper is organized as follows: In Sec. II we intro- 
duce the Euclidean action to describe the quantum decay 
in DJTL's with dissipation. In Sec. Ill, we discuss the 
saddle-point solutions of the action as a function of the 
temperature and the driving current. An iterative per- 
turbation procedure to calculate these extrema close to 
a saddle-point bifurcation is presented in Sec. IV. This 
method is then applied to evaluate the split-instanton so- 



lutions in Sec. V. The various decay regimes, which are 
summarized in a diagram and the corresponding relax- 
ation rates are discussed in Sec. VI. Finally, the conclu- 
sions are drawn in Sec. VII. 



II. MODEL 

A. Decay rate 

The probability per unit time for a metastable state 
to decay at a finite temperature T is related to 
the imaginary part of the free energy by F = 
{2/h)lmF, as was shown by Affleck. □ In the path- 
integral representation, the free energy is given by 
F = -kT\n[§V[q] exp(-S'B[q]/?i)], where Se is the Eu- 
clidean action and q is a vector representing the DOF. 
If the energy barrier that the system has to overcome in 
order to leave the metastable state is large, the sum over 
paths is dominated by the extremal trajectories and the 
latter can be treated in the semiclassical approximation. 
The decay rate then reads 

r = Ae-^/'', (1) 

where B is the Euclidean action Se evaluated at the ex- 
tremal saddle-point trajectory and A is a prefactor that 
can be obtained by considering the fluctuations around 
the extremal path. In this work, we will determine the 
exponent B. 



B. Euclidean action 

The Euclidean action of a system of A^ identical JJ's 
in the presence of a bias current / is 

Se^ dt [K{^) + V{v) + D{<f)] , 
Jo 

where t is the imaginary time, (p = {(po, . . . , tfN-i) repre- 
sents the phase differences across the junctions, V is the 
potential energy, and D models the dissipation. The "ki- 
netic" energy containing the capacitances C of the JJ's 
is given by 




where m = C(4'o/27r)^ corresponds to the "mass" of a 
fictitious particle and $0 = hc/2e is the flux quantum. 
The potential energy consists of two parts, 

where U{(fi) — ^[1 — cos(iy9„) — (//A^/c)(p„] represents 
the tilted washboard potential of the driven JJ's that 
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arises due to the relation between currents and gauge in- 
variant phases across the junctions. Here / is the total 
current through the system, Ic is the critical current of 
a single junction, and Ej — {^o/2tt)Ic is the Josephson 
energy. We concentrate on the experimentally most inter- 
esting limit of currents close to criticality, NIc — I <^ Ic, 
where the tilted washboard potential can be well approx- 
imated by its cubic expansion, 



N-l 



N-l 



n=0 



Here e — — I /NIc) is a small parameter that indi- 

cates the distance from criticality / = NIc, 'f = 7r/2 — eis 
the value of the phase at the minimum of the potential, 
and Uq is an irrelevant constant that will be dropped 
in the following. Taking only the self-inductances L 
of the loons-, into account and neglecting the mutual 
inductancesEil, the interaction energy between the loops 



A/'(q) = -1790 - '2<loJ2 ~ ^ J2 {<l2k ~ q2{N-k)) 



k=l 



N-l 



2^ ^ ^ ^ 

g— 2^ QmQk [Qm+k + Qm-k — Q2N-m-k) 

ni>k—l 

with gfc = for fc > - 1. 

The action in terms of the dimensionless coordinates 
reads 



Se ^ g I drL 



(5) 



where the Euclidean Lagrangian is 



Le = V(q) - 



e dq 



ny/j dr J _^ dr' 



Eif ) = ^ ^ ((^„+i - ^^„)^ 

where /3 = LI^/Ej is the McCumber parameter. We 
model dissipation due to an OhmiCpShunt resistance R 
using the Caldeira-Leggett approachO, 
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Here g — NUEje^ /ttT 1 is the semiclassical param- 
eter, 9 — ■Kr](3T/2hEjsm'^{n/2N) is the dimensionless 
temperature, and J = (/3e)^/[2 sin(7r/2A^)]'* is the di- 
mensionless current. Note that J = corresponds to 

/ ^ NIr. 



^^a^'a^^T^" 



III. EXTREMA OF THE EUCLIDEAN ACTION 



where rj = l/R is the phenomenological friction coef- 
ficient. In the following we will treat the overdamped 
limit and neglect the contribution of the capacitive term 
K{ip) to the action. It is convenient to perform a trans- 
formation to dimensionless normal coordinates. 



Vn{t) = + 2ego(T) 



-2V2 e ^ qk{r) cos 



fc=i 



■nk{n+ 1/2) 
N 



(2) 



Note that the qk are functions of the dimensionless imag- 
inary time r = ^TxTtjfi — tt. We now define the dimen- 
sionless potential energy. 
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where 



fc=0 



9 / Trfc , 
2 \ \ 



LLu = — sm' , 
^ I3e V 2N 



(4) 



are the eigenvalues of (9ma„V) evaluated at the local 
minimum q = 0. A/'(q) contains the cubic terms 



By applying the variational principle to the Euclidean 
action, one finds the classical equations of motion in 
imaginary time 



V,V(q) 



.1 + 
dr — — cot 

or 



0. (6) 



Their solutions are given by the extremal trajectories, 
of which the saddle-point solutions are of special inter- 
est, since they lead to the decay of the chain from its 
metastable state. At high temperatures, quantum fluctu- 
ations play a minor role and the solutions of the equation 
of motion are time independent, dr<\ = 0. Hence, they 
are given by the extrema of the potential, VqV(q) = 0. 
However, below a crossover temperature 9q, quantum 
tunneling becomes relevant for the decay process and 
the solutions of the equation of motion are a function 
of the imaginary time. In the following paragraphs, we 
first analyze the extrema of the potential, which deter- 
mine the decay in the thermal regime. Then we derive 
an expression for the crossover temperature from ther- 
mal to quantum tunneling and finally analyze the time- 
dependent extremal solutions of the Euclidean action, the 
instantons which lead to decay via quantum tunneling. 
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A. Saddle points of the potential energy 

A trivial type of saddle-point solution can be readily 
constructed from physical arguments. Consider the case 
where the attractive interaction between the particles is 
much larger than the energy barrier. At high tempera- 
tures, the strongly coupled particles most probably are 
thermally activated over the barrier all at once and the 
chain basically behaves as a rigid rod. In this case the 
saddle point of the potential V is identical to the lo- 
cal maximum of the single DOF potential, (po — . . . — 
ipN-i = e + 7''/2, which reads qq — 1, qk>o = in normal 
mode representation. If, on the other hand, the energy 
barrier is of the order of the interaction strength or even 
larger, another saddle point of the potential emerges. 
Then above a certain barrier height the chain preferably 
decays via a kinked saddle-point solution with qk>o 7^ 
that we call elastic. 

Taking the second derivatives of Eq. (j^) at q^s = 
(1,0,..., 0) we evaluate the eigenvalues of the curvature 
matrix at the thermal rigid saddle point. 



Ar=a,^V(q.,) = -sin^( J^)-2 



I3e 



2N 



(7) 



Note that A^* — 2(1/VJ — 1) becomes negative for J > 1, 
indicating that another saddle appears. This elastic so- 
lution has a lower activation energy and hence becomes 
the most probable configuration leading to escape from 
the metastable stateE 

To determine the saddle points, one has to solve a sys- 
tem of N coupled nonlinear equations. So far, elastic 
saddles have been calculated exactly for N = 2, Ref. [6] 
and = 3, Ref. [9]. Close to the crossover J>1, the 
Qk>o are small. Hence, for A^ > 3, approximate solu- 
tions can be found.l3 Expanding around the rigid saddle, 
go = 1 + 9o and qk — qk for fc > 0, we approximate the 
potential energy by 



V 



k=0 



rs~2 
k Ik 
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(2go 



\f2> 
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(8) 



Solving VgV = 0, one finds ql = A5''A5;''A^V(4A5''-h8A^^), 
qo = 2ql/\\f, q2 = \/2q{l\Tf , and qk>2 = 0. Hence the 
eigenvalues of the curvature matrix evaluated at the elas- 



tic saddle for /c ^ 1 are A^ 



AI 



whereas Af = 2|A5;^|. 
Inserting these results into Eq. (||), one finds that the 
activation energy is reduced due to the elasticity. Close 
to the crossover, where J>1, it is given by 



v4 



4(a;s« + 2A-) 



1 

71 



(9) 



B. Crossover temperature from thermal to quantum 

decay 

In the following, we derive the crossover temperature 
^0- Since the time-dependent quantum saddle-point so- 
lutions are periodic in r with a period 27r, they can 
be represented by a Fourier series. Near the crossover, 
the Fourier series can be well approximated by q(T) w 
qts + pcos(r), where qts is the thermal saddle point un- 
der consideration (rigid or elastic) and p is a small cor- 
rection term due to quantum fluctuations. Substituting 
q(T) into the linearized equation of motion. 

Op 



[dkdiV (qts)] PCOs{t) 



dr' sin(r') cot 



one finds an eigenvalue equation for the potential curva- 
ture matrix, 



[dkdiV{qts)] P 



29 

71' 



The only negative eigenvalue of the curvature matrix 
[dkdiV{qts)] evaluated at the saddle point of the potential 
is given by Aq, the curvature along the unstable direction. 
Hence the crossover temperature is given by 



|Ao|\/j 



(10) 



For the rigid regime, where X^^ = —2 we thus find 

0o(J<l) = \/j, (11) 
and in the elastic regime with Xff « A" — 2A"A2*/(Aq* — 



1 



2cos2(7r/27V) 



(12) 



C. Instantons 

For 9 < 60 quantum tunneling becomes relevant and 
the instanton solutions dominate the decay from metasta- 
bility. For J > 1, the crossover from thermal to quantum 
decay is of second order. A detailed procedure show- 
ing how to obtain the saddle-point solutions close to the 
thermal-to-quantum crossover was given in Ref. [12] and 
can be applied to calculate the quantum elastic solutions. 
In this paper we will concentrate on the current regime 
J < 1, where in addition to the transition from the ther- 
mal to the quantum rigid regime a crossover from the 
quantum rigid to the quantum elastic phase can take 
place. 

The rigid quantum solution is found by setting qk — 
for A: > 0. In this case the equations SSs/Sqk — with 
k > are trivially satisfied and the remaining equa- 
tion for k — describes the thermally assisted quan- 
tum tunneling of a single degree of freedom go • Its solu- 
tion is the W|ell-known instanton obtained by Larkin and 
OvchinnikovQ 
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0/ 1 - ^1 -6l2/6'o cos(r) 



(13) 



The operators Qk are defined as 



Inserting q^'' and q^^y^ — into Eq. (||), one obtains the 
extremal action in the rigid quantum regime 



Bqr — Bq 



^ 3 Uo 



(14) 



where Bq — 2TrNrie'^. 

As for J > 1, nonuniform saddle-point solutions of the 
action exist in the quantum regime in a certain param- 
eter range even for J < 1. If the action evaluated at 
this extremum is lower than Bqr, the nonuniform config- 
uration is the most probable one leading to decay from 
the metastable state via quantum tunneling. Tuning 
the temperature at a fixed bias current, a nonuniform 
saddle-point solution can develop in two different ways. 
One possibility is that a less probable nonuniform config- 
uration, which coexists with the rigid saddle point above 
a critical temperature 0i , becomes the lowest-lying saddle 
point of the Euclidean action below 9i. Then the most 
probable configuration abruptly changes from uniform to 
nonuniform. Since the first derivative of the rate dgT is 
discontinuous at 9i, the crossover is of first order. An- 
other scenario is encountered, if at a critical temperature 
02 the rigid saddle point bifurcates into new saddle points 
which have the lowfgt action. This crossover, known as 
instanton splitting,Lli3 is of second order. 

A strategy to determine nonuniform saddle-point so- 
lutions for J < 1 is to first search for a saddle-point 
splitting and then to verify whether a first-order transi- 
tion might have occurred before the bifurcation has taken 
place. If a first-order transition can be ruled out, the new 
bifurcated saddle points have the lowest action. In this 
case the bifurcation causes a second-order crossover from 
a single to a split-instanton regime. 

Following this idea, we first identify the saddle- 
point bifurcation, calculate the split instantons and test 
whether or not a first-order transition has already oc- 
curred. 



IV. ITERATIVE PERTURBATION SCHEME 

In this section we present an iteration scheme to calcu- 
late the split-instanton solutions for J < 1. We start by 
expanding the coordinates around the single instanton 
solution, 



(15) 



and rewrite the Euclidean action in terms of the new 
variables, 



r 1 

Se ^ Bqr + g d,T qkQkqk + A/'(q) 

•^-^ L fc=o 



(16) 



Qkqkir) = - 

9 



oqk(T)dqk[T') 



dr' 



To determine the split instantons, we have to find saddle 
points of the action with nonzero qk- Hence we have to 
solve the equations of motion. 



Qkqk 



9AA(q) 
dqk 



(17) 



which constitute a system of coupled nonlinear differ- 
ential equations. In general, they cannot be solved ex- 
actly. However, close to the saddle-point bifurcation, the 
extremal amplitudes qk are small and we can calculate 
approximate solutions by applying an iterative perturba- 
tion scheme. This leads to a hierarchy of inhomogeneous 
linear equations. 



Qkqt^ 



F, 



ii) 



(18) 



where i denotes the iteration step. In the first iteration 
i — \ we take only terms into account that are linear 
in qk- The higher-order terms on the right-hand side of 
Eq. (0) are neglected. Thus, F^^^ = and Eqs. (|l|), 
which have to be solved, are homogeneous. For z > 1, 
the amplitudes calculated in the previous iteration are 
substituted into dN such that the inhomogeneous terms 
are given by 



9AA[q('-i)(r)] 
dqk 



(19) 



After each iteration step i, we thus obtain approximate 
(special) solutions for the amplitudes qk by formally in- 
verting Eq. (Iisl), 



(20) 



Of course, a straightforward inversion is not possible, if 
Qk is singular. Below, we will discuss how to handle 
equations with a singular operator. 

The inversion is most conveniently performed by rep- 
resenting Eq. (|o|) in terms of the eigenfunctions of the 
operators Qk which we will determine now. One realizes 
that the operators Qo and Qk only differ by a constant 
term 

Qk = Qo + — Mo- 

They trivially commute and have a common set of eigen- 
functions ilim- The eigenvalues Af„ of the operators Qk 
are related by 
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A* 



A„, + /ife - Mo- 



lt is, therefore, sufRcient to concentrate on the eigenvalue 
problem 

which was studied by Larkin and Ovchinnikovil in 
the context of single-particle tunneling with dissipation. 
They obtained the spectrum 

A? = -2 + 2ae, 
K = -2ac, 
A^ = 0, 

AO,-2[l + (|m|-2)0/0o], |m|>2, 



where ac = 1/2 + a/5/4 - 9'^/9l and showed that the 
eigenfunctions 



n— — oo 



have Fourier coefficients of the form 

Bm{Cm,n + C?m,n)j 'T- ^ 0, 



m,n+dm,«), ±m > 0, n < 0, 



with dm,n = for \m\ < 2 and |n| + 2 > |m| > 2. Note 
that the ipm are even (odd) for positive (negative) m and 
the Brn are chosen such that the eigenfunctions are nor- 
malized, 



(V'm,V'm) = / dr ipl^ir) = 1. 

J — TT 

For m = 0, 1 they obtained 

C„.,„ - (|n| - ^Al^ e-^l"! (21) 



with tanh5 = 9/9o. In fact, calculating the remaining 
coefficients for m < — 1 and m > 2, we find that Eq. ( |2l| ) 
holds for any m with 



\n\ < \m\-2, 
71 = m + 2 = 0, 

n = m — 2 — 0, 



i(S-l)e-''('"l-4), |n| = |m|-2>0. 



With these results, we now represent Eq. ( po| ) in terms 
of the basis tpm with 



m— — oo 

and obtain a special solution in terms of the coefficients 

41 = /il/Af„, if At / 0. (22) 

If for some k' and m' the eigenvalue A^, = 0, the op- 
erator Qk is singular and a unique solution of Eq. ( [T^ ) 
cannot be found within the «th iteration. However, after 
performing all necessary iterations, the solutions have to 
be the lowest-lying saddle point of the Euclidean action. 
This constraint enables us to determine the so far arbi- 
trary coefficients c)./^, by requiring that the Euclidean 
action as a function of the coefficients has to be minimal. 
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V. NONUNIFORM INSTANTON SOLUTION 

After we have explained in detail how to obtain the 
approximate solutions, we are now ready to perform 
the calculations explicitly. Let us define the parameter 
a = (/ii — /io)/2 = l/'s/j. First, we show that the in- 
stanton splitting occurs at a = Ofc and then we apply the 
perturbation scheme to determine the split-instanton so- 
lutions. 

The negative Ag indicates that the operator Qq has 
an unstable mode, which is responsible for the imaginary 
part of the free energy and hence for the finite decay rate 
of the metastable state. For a > the spectrum of Qi 
is positive definite. Since /Xfe — /iq > 2a, for fc > 1 A^ > 
and hence all higher modes qk>o are stable. For a < ac, 
the lowest eigenvalue Aj of Qi becomes negative, indicat- 
ing that the corresponding mode also becomes unstable 
and that a new saddle point with a lower Se exists. Thus 
at a = ac a split-instanton solution emerges.u 

To determine the split-instanton solutions for a<ac, 
we now apply the iterative procedure and solve Eqs. (|lq). 
In the first iteration Fm ^ — 0. According to Eq. (|2 



most of the coefficients 



.(1) 

''k^m 



are zero except c^q'Li and 



c^^Q. The coefficient Cq^^^ cannot be uniquely determined 
since A'^i = 0. However, the corresponding odd eigen- 
function ip^i is associated to imaginary time translation 
symmetry and does not contribute to the value of Se- 



We have, therefore, the freedom to choose c, 



0,-1 







m— — oo 



within this and the following iterations. At a = ac, 
where Aq = 0, the operator Qi becomes singular. Here 

the instanton splits since the coefficient ^ = c^^g of the 

dangerous mode -00 of Qi can have a finite value that re- 
mains to be determined by minimizing SEiC)- To lowest 
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order, the split-instanton solution at a 
by 



OLc is thus given 



0. 



(23) 



In analogy with the Landau theory of phase transitions 
we can interprete Se as a thermodynamic potential and ^ 
as an order parameter. A finite ^ indicates the existence 
of a quantum elastic solution. 

Using Qi{a) = Qi(q^c) + 2{a — ac), recalhng that 
Qi(ac)^i = and inserting the perturbative result (|2^ ) 



into Eq. (16), one obtains the split-instanton action up 
to terms quadratic in the dangerous mode, Se{C) = 
Bqr + g{a — ac)C^- However, in order to be able to min- 



imize the action as a function of C, for a < ac 
the terms quartic in C have to be calculated, 

SE{0 = B,r+g[{a-ac)e+5C^]. (24) 

The case 5 < 0, which indicates that a first-order tran- 
sition has occurred will be discussed later in more de- 
tail. For (5 > 0, the minimal value of Se is given by 
Id = {etc — a)/2(5 and the extremized action reads 



Bqe — Se{0 — Bq 



4(5 



(25) 



Note that C, is small close to ac and can be regarded as 
a perturbation parameter. 

In the first iteration we considered corrections to the 
single instanton solution of the order C,. In order to de- 
termine J, the split-instanton solution up to orders of 



has to be treated. Consequently, we have to per- 
form the second iteration of the perturbation procedure. 

The inhomogeneous terms F^^ in Eq. (|l^) are found 



by substituting q^. 



F, 



(2) _ 



0, hence q'^^ = ql}^ 



into Eq. (|l9|). For k ^ 0,2 the 
whereas for /c — 0, 2 we 



obtain F^"' = 2C^{tPo)'^ and F^'> = V2C^(V'o)^- The 
remaining task within this iteration is to solve the equa- 
tions 



QoqP = 2C2(Vo)' 
Q2gf ) = V2C^{^o 



at least Representing (ipo) in the basis t/jm, we obtain 



E 



Note that the odd ij^m with m < do not appear in the 
sum, since (V'o)^ is an even function. The coefficients 
are given by 



Ln 



For our purposes it will be more than sufhcient to con- 
sider only the first three coefficients ao,ai and 02, since 
o-m/o-m+i ^ 1 in the entire quantum regime. After te- 
dious but straightforward calculations, we find 



ao 



ai 



02 



' 13, 



/27r 



/2Tr 



17\ 

4 } 



'37, 



19^ 
2 J 



'69 



8^ 



Oic 



\ 21 
) 9^ 



{2ac 



3\ "5 
2) W 



{ac 



So 



3/2 



00 



63 



-2ac 



tOIc 



7\ «o 
2) W 

13 

8 



1/2 r 



(2ac + i) ^ - («c + i) 



So 



J e 



tOIc 



) W 



1 f^o 
4 P 



2W 



1/2 



m > 2. 



(2«c + I) ^ - («c + I) 



(26) 

(27) 

(28) 
(29) 



Note that for 9 ^ 0q the coefficients converge to 
ao 1/^/2tt and ^ for fc > 0. Substituting 
fo,m = 2C^a„ and f2.m = \/2C^a,„ into Eq. (||) we 
find excellent approximations of the expansion coeffi- 
cients Cg^^ and c^2m- Having determined qj^"^ and q^'' , the 
second iteration of the perturbation procedure is com- 
pleted. 



I 

Evaluating Eq. ( p6| ) with q-Q^'' , q^^"^ oc i we obtain the 
coefficient S of the quartic term in S'e(C), 

lyrlMr)? I^qi'^ + (30) 
Performing the integral and using the orthogonality of 
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FIG. 1. The fourth-order expansion coefficient S of the 
Euclidean action Se{0 = Bqr + g [{a - ac)C + + 7C'^] 
as a function of dimensionless temperature 9/9q for A'' = 2, 3 
degrees of freedom and the continuous hmit N oo. 
In the low-temperature limit 5 is positive and diverges as 
S{9) ^1/9. At the tricritical temperature 9tc it vanishes, 
5{9tc) = and becomes negative for 9 > 9tc- The negative 5 
indicates a first-order transition. 



the Tpm, we find 



m=0 



2 



1 

A^" 



(31) 



The function (5(6'/6'o) is shown in Fig. for A^ = 2, 3 
and N —f oo. At first it may be surprising that the 
curve for N ~^ oo lies in between the curves for A^ = 2 
and N = 3. The reason is that for N — 2 the 
mode q2 does not exist and S = — X]^m/^m- Since 
A^ = AJ'^ + 8acCOs'^{TT/2N) > for A^ > 2 it is clear 
that 6{e/0o,N = 2) > S{0/eo,N > 2). When eval- 
uating Eq. jsTl ) for > 3 one obtains the relation 
S{9/6o,N) < S{e/eo,N+l). In other words, for iV > 3 
the graphs lie in between the ones for A'^ = 3 and N oo. 
In the limit 9 we find d oc Oq/O when taking into ac- 
count only the leading terms in Eqs. (p6|)-(p8|). Hence, 
g/6 in Eq. ( |25| ) converges to a constant value for 9^0. 
With increasing 9 the coefficient 6 decreases and vanishes 
at the characteristic temperature 9tc, where Eq. ( p5| ) 
looses its validity. At 9tc the first-order and second-order 
transition lines merge. In analogy to the classical theor; 
of phase transitions it is called tricritical temperatureBI 
Above 9tc the parameter 6 becomes negative, indicating 
that the transition from rigid quantum to elastic quan- 
tum decay becomes first-order-like. The values of 9tc are 
given in Table | for A^ = 2, 3 and N ^ oo. Note that 9tc 
is smallest for A^ = 3 and increases monotonically with 
AT for A^ > 3. Recall that 9tc is largest for A^ = 2 due to 
the absence of q2 in a system with only two DOF. 

We now concentrate on the case where i5 < 0. Then 
in order to find the minimal values of Se, we have to 
determine terms of the action cx , 



/V 

i V 


9 


■3 
O 


oo 


Jtc 


0.8424 


0.7652 


0.7938 


9tc 


0.8719 


0.8000 


0.8275 


S'{9tJ9o) 


-2.274 


-1.384 


-1.675 


l{9tc/9o) 


0.4118 


0.2105 


0.2698 



TABLE I. Numerical values of the tricritical current Jtc , 
tricritical temperature 9tc, the derivative of the coefficient of 
the forth-order term S'{9tc/9o) and the coefficient ^{9tc/9o) 
of the sixth-order term for various numbers of degrees of 
freedom TV. 



which are obtained in the third iteration (i = 3) of the 

perturbation procedure. For A; ^ 1, 3 one has q'j^^ = q^j^^ ■ 
Inverting the equations 



9f ) CMr), 



QigP^ = (4g^'^ + 2V2, 
Oagf ^ = 2V2C^o(t)#^ 

numerically and inserting the values into 



1 



(3) 



V2 



{-V2C^o(r)(zf^g| 
(Mr)] [2q^o'^ + V2qi'^] (Mr)}, 



(33) 



we calculate the coefficient 7(0/6*0). Minimizing Eq. (|3£ 
with respect to ( for a — ac < (5^/37, one obtains, in 
addition to the rigid instanton solution with ^ = 0, an 
elastic instanton solution with 




1 - ^(a - ac) 



(34) 



Se{0 = B,r +g[{a- a,)^ + SC + iC] 



(32) 



The first-order transition occurs at a = ai = ac + 5'^ /Aj 
when the nonuniform solution becomes the global mini- 
mum of the action, SsiC) = Se{0) — 0. 

Using the perturbation scheme, one could, in principle, 
determine the split-instanton solution to arbitrary order 
in C. For our discussion of the behavior of the decay 
rate close to the crossover from the single instanton to 
the split-instanton regime, the calculation shown above 
is sufficient. 



VI. RESULTS AND DISCUSSIONS 



In this section, we discuss the various decay regimes 
which are presented in the decay diagram (Fig. ||). Let 
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FIG. 2. The decay diagram for iV = 3 de- 
grees of freedom in terms of the dimensionless current 
J = P^{l-I/NIc)/[8 sin'*(7r/Af)] and the dimensionless tem- 
perature 9 = TvriPT/2hEjsm^{-K/2N). Criticahty I = NIc 
corresponds to J = 0. Above 9o{J) the system is in the 
thermal regime, and at J = 1 a second-order crossover from 
rigid to elastic thermal decay occurs. At 6*0 (J) a second-order 
crossover from thermal to quantum decay takes place. The 
quantum regime 8 < 9o is again separated into rigid and 
elastic decay. For Jo < J < Jtc the crossover from rigid to 
elastic decay at 92{J) is of second order. For Jtc < J < 1 
the crossover indicated by 9i{J) is first-order like. Though 
the diagram is similar for different A*', the temperatures 9tc, 
9i{J), a,nd9o{J > 1) are altered. For comparison, the dashed 
curve shows 9q{J > 1) for a dc SQUID (A = 2). 



US start with the thermal regime 9 > Oo{J), where the 
decay occurs via thermal activation. For J < 1, the 
coupled DOF behave like a single DOF since the cou- 
pling energy is large compared to the thermal or the 
barrier energy. Then the system is in the rigid thermal 
regime^ Increasing J, one enters the elastic thermal decay 
rcgimeo passing the second-order crossover line at J = 1. 
On the other hand, starting in the thermal rigid phase 
and reducing 9, quantum fluctuations become important 
and at Oq{J) a second-orde|r crossover from thermal to 
quantum decay takes place.a Two characteristic currents, 
Jo = 0.3820 and Jtc as given in Table | become important 
in the quantum regime. Below Jq, the system preferably 
decays via the single instanton or rigid quantum saddle- 
point solution. For J > Jq a transition from the rigid 
to the elastic quantum decay regime becomes possible. 
For Jo < J < Jtc, the crossover is of second order and 
is caused by a saddle-point bifurcation of the Euclidean 
action occurring at a = ac. The dimensionless crossover 
temperature is then given by 



= fj + VJ- 1 



1/2 



(35) 



For Jfc < J < 1, the crossover is of first order. The 
transition occurs at a = ai. Near 9tc, we approximate 
S « S'{9tc/9o)[{9 - 9tc)/9o] and find that in this limit the 
crossover line is given by 



(36) 



The numerical values for 7, 9tc, and S'{9tc/9Q) are given 
in Table |. Note that since da9i diverges as a ac, the 
slope of 9i{J) is infinite at the tricritical point. For J > 1 
the transition from the thermal to the quantum elastic 
region is again of second order. The crossover tempera- 
ture is then given by 9o{J > 1) [see Eq. (p^J. The decay 
diagram of an overdamped JJ array with N = 3 junc- 
tions presented here is similar to that of the dc SQUID 
with N — 2. Qualitatively, the diagrams exhibit the 
same features for all A^. The reason is that for J < 1, 
the transition lines 9q{J) and 92 (J) are determined by 

the long- wavelength modes ^q'^'' and qi , respectively, and 
hence are independent of A^. However, there is a differ- 
ence between the diagrams on a quantitative level, since 
9tc, 9i{J) and 9o{J > 1) are parametrized by A^. For ex- 
ample, compared to the dc SQUID, the first-order tran- 
sition region is enlarged for N > 2 and is the largest for 
N = 3. 

The remaining task is to discuss the decay rate F ~ 
exp(— _B/?i) in the four regimes. To exponential accuracy, 
F is determined by the extremal of the action B, which 
is given by the Euclidean action Se evaluated at the rel- 
evant saddle-point configuration (f^, B{J,9) — Se[^s\- 
The behavior of the rate in the thermal regime was 
discussed in Ref. [9]. Since the thermal saddle points 
iptg are independent of imaginary time, B — hV{ipg)/T 
and the rate reduces to the classical Arrhenius form, 
F ^ exp{—V{<f^)/T). In the thermal rigid regime B 
is given by 



B 



2Bo y/j 



tr 



J < 1. 



(37) 



Realizing that in the thermal regime S = — (2/Aq* + 
l/\^')/{2n) and recalling that a = l/Vl, one finds, with 
the help of Eq. (H), the thermal elastic result. 



Bte ^ Btr - 9 ^"^ J>1. 



(38) 



In the rigid quantum regime, the action Bqr is given by 
Eq. (|l|). Inserting Eq. (H) into Eq. (|3|) we find the 
extremal action in the quantum elastic regime for J < 1, 



a — OLr 



2772 

+2 [<52_37(a-ae)]'/'} 



26-' 



(39) 



where 5{9/9q) is given by Eq. (|3ll ) and ^{9/9q) was cal- 
culated numerically using Eq. (|33) . The rates for various 
currents J are displayed in Fig. H as a function of tem- 
perature. For J < Jo the system is in the rigid regime 
for all temperatures 9 (see Fig. ^). For 9/9q > 1 the 
thermal rigid result (p^) applies. In the rigid quantum 
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FIG. 3. The extremal of the action B as a function of 
temperature 6 for various normalized currents J. The dashed 
hne shows the purely thermal behavior in the rigid regime. 
The rigid quantum result is represented by the solid line. 
The dotted and the dashed-dotted lines display _B of a sys- 
tem with N = 3 degrees of freedom for J = 0.6 and J — 0.9, 
respectively. For J = 0.6 a second-order crossover to the 
split-instanton regime occurs at 62 ~ 0.790^0. A first-order 
crossover takes place for J = 0.9 at 61 = 0. 975^0. The inset 
shows the cusplike shape of B close to 9i. 

regime O/Oq < 1, in comparison with the purely thermal 
result, the rate is increased due to quantum fluctuations 
according to Eq. (p|). In the chosen representation Bqr 
is independent of system-specific parameters. Experi- 
mentally measured decay rates of rigid systems should 
thus collapse onto one curve. For J > Jq tunneling of 
nonuniform instantons becomes possible and B = Bqg is 
reduced further compared to Bqr- In Fig. |^ we displayed 
B for a system with = 3 degrees of freedom. As an 
example for the behavior of the rate close to a second- 
order crossover to the split-instanton regime we chose 
J = 0.6 < Jtc- The crossover occurs at 6*2 = 0.7906'o. For 
J = 0.9 > Jtc, the behavior of the rate is different. The 
slope of Bqe changes abruptly at 9i = 0.975^0, which 
indicates the occurrence of a first-order crossover. 



VII. CONCLUSIONS 

In the present work, we studied the decay of 
metastable states in current-driven parallel coupled one- 
dimensionsal Josephson-junction arrays at zero voltage 
in the overdamped limit. We model this system by N 
elastically coupled DOF trapped in the minimum of the 
single-particle potential and interacting with a bath of 
harmonic oscillators. The escape from the trap can be 
induced by thermal or quantum fluctuations. Three en- 
ergy scales determine the decay behavior of the system; 
the temperature, the barrier height of the trap and the 
interaction between the particles. Accordingly, one finds 
four different regimes for the decay rate which we summa- 
rized in a decay diagram in Fig. |[ To calculate the decay 



rate we use the thermodynamic method. In the saddle- 
point approximation, the decay is then determined by 
the most probable configurations leading to an escape 
from the trap which are given by the saddle points of the 
Euclidean action. In the thermal regime the saddle-point 
solutions are independent of the imaginary time and iden- 
tical to the saddle points of the potential energy.B If the 
interaction between the DOF is strong compared to the 
barrier energy, rigid configurations dominate the decay. 
Reducing the bias current, the barrier becomes larger and 
above a critical value the system preferably decays via an 
elastic configuration. On the other hand, starting in the 
rigid thermal regime and lowering the temperature, quan- 
tum fluctuations become important and the decay most 
probably occurs via the rigid quantum saddle-point or 
single-instanton solution of the Euclidean action.El Inside 
the quantum region, an elastic regime can again be en- 
tered by increasing the barrier above a critical value. In 
order to determine the nonuniform instanton solutions 
of the quantum elastic regime, we worked out an itera- 
tive perturbation procedure. We performed the calcula- 
tions close to the crossover from rigid to elastic quantum 
decay analytically up to second order and realized the 
third-order calculation numerically. We were then able 
to give quantitative results for the decay rate including 
the quantum elastic regime. The behavior of the decay 
rate is similar for SQUID's, DJTL's, and long JJ's. In 
the rigid regime the decay occurs via a saddle point that 
is uniform in space and hence, the qualitative nature of 
thermal or quantum decay is not sensitive to the number 
of DOF. Further, the crossover from rigid to elastic de- 
cay is caused by the excitation of long-wavelength normal 
modes of the system, which are equivalent in the three 
physical systems discussed here. 

We want to emphasize that although our conclusions 
are drawn for overdamped systems, the reasoning and the 
procedure also apply for the underdamped case. Indeed, 
on a qualitative level, the understanding of the quan- 
tum rigid-to-elastic crossover in underdamped DJTL's 
and long JJ's can be obtained jjn the basis of the theo- 
retical work on SQUID's (Ref.l3). However, in order to 
have quantitative results, one has to extend the theory 
following the scheme proposed here. 

One interesting aspect of the quantum rigid-to-elastic 
crossover is that depending on the current, it can be ei- 
ther of first or second order, whereas all other crossovers 
that we discussed are of second order. Even more fasci- 
nating is the fact that this crossover is an intrinsic quan- 
tum property and can be regarded as one further evi- 
dence for MQT, if measured. Experimental verifications 
of the predicted enhancement of the decay rate due to 
the elastic properties (see Fig. ||) are thus highly desir- 
able. An experimental detection of the first-order-like 
crossover would be challenging, but seems to be difficult, 
because the cusplike behavior of the rate at the crossover 
is not very pronounced and occurs in a small current 
interval. In standard experiments, the rate is obtained 
from the switching current histogram. The current in- 
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tervals of the histogram have to be much smaller than 
8iV/c(l - Jtc)sin''(7r/2iV)//32 and the number of events 
per interval large in order to resolve the cusplike feature. 
Furthermore, it would be convenient to perform the mea- 
surements on systems with = 3 DOF since in this case 
the first-order region is the largest (see Fig. ||) . 

In sum, we calculated the decay rate of overdamped 
current-biased one-dimensional Josephson-junction ar- 
rays at zero voltage (including SQUID's, DJTL's, and 
long JJ's) analytically and numerically in several distinct 
decay regimes. An experimental observation of the pre- 
dicted enhancement of the decay rate in the elastic quan- 
tum regime would give further evidence for macroscopic 
quantum tunneling in these systems. 
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